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The geometrical mapping of algebraic nuclear cluster models is investigated within the coherent 
state formalism. Two models are considered: the Semimicroscopic Algebraic Cluster Model (SACM) 
and the Phenomenological Algebraic Cluster Model (PACM), which is a special limit of the SACM. 
The SACM strictly observes the Pauli exclusion principle while the PACM does not. The discussion 
of the SACM is adapted to the coherent state formalism by introducing the new 50(3) dynamical 
symmetry limit and third-order interaction terms in the Hamiltonian. The potential energy surface 
is constructed in both models and it is found that the effects of the Pauli principle can be simulated 
by higher-order interaction terms in the PACM. The present study is also meant to serve as a 
starting point for investigating phase transitions in the two algebraic cluster models. 
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I. INTRODUCTION 

Phase transitions in algebraic models have been dis- 
cussed since the 1970's 0J. The principal method is the 
use of coherent states [2j • The expectation value of the 
algebraic Hamiltonian, with respect to a coherent state, 
is defined as the semi-classical potential. Through the be- 
havior of that potential, as a function in the parameter 
space, phase transitions and their order can be studied. 
The basic description, applied to the IBA Q, was pre- 
sented in 0] . The method of coherent states also delivers 
a geometrical mapping of the algebraic model in consider- 
ation, providing in this manner an easy interpretation of 
the dynamical symmetry limits. Other studies on phase 
transitions in the IBA are published in Q. 

In 0], a complete classification of phase transitions in 
algebraic models is presented, restricted to Hamiltonians 
with up to two-body interactions. For the vibron models 
^,0] the transitions turned out to be of second order. In 
8| , a second order phase transition was also encountered. 
In the U (3) boson model was studied and it was noted 
that for very large number of bosons, the transition may 
turn over into one of first order. The transitions were 
investigated using the overlap of the ground state with 
that of the 0(2) limit and searching for a step like behav- 
ior. In other words, no discontinuities of the derivatives 
of the potential were considered. 

Coherent states have not only been applied to the IBA 
or atomic molecules @, 0], but also to other algebraic 
models, which have a microscopic origin within the shell 
model. In [ic| the geometrical mapping, using the vec- 
tor coherent state method I 111 was applied, mapping the 
pseudo-sym plec tic model [12[ to the geometric model of 
the nucleus [13j . The geometrical mapping turned out to 
be very useful in calculations of nuclear spectra 14. Il5| 
and predicting the spectra of super heavy nuclei 16[ . In 



[1 71 ] , the coherent state method was used to obtain a geo- 
metrical mapping of the Semimicroscopic Algebraic Clus- 
ter model (SACM) [H,[I1]. The SACM is an algebraic 
cluster model which takes into account the Pauli exclu- 
sion principle. It allows for investigation of the effects of 
the Pauli exclusion principle on the geometric potentials 
and on the order of the phase transition. Phase transi- 
tions in single nuclear systems were also investigated in 
0, HHi related to the symplectic model of the nucleus. 

Note that phase transitions in nuclei, though not ex- 
plicitly stated as such, were already studied in 1972 in 
the first edition of [l3| . A standard curve discussion was 
applied, while a phase transition was denoted as a shape 
transition. In the recent treatments, the main difference 
is the classification in terms of the order of phase transi- 
tions. 

Since its initial use, interest in the geometrical map- 
ping, using the coherent state method, has not been lost. 
One of the main reasons is that multiple particle sys- 
tems in conjunction with phase transitions can be easily 
treated. These methods can provide important insight 
into how to treat, in general, complicated many-body 
systems. 

In recent years, phase transitions in atomic and nu- 
clear molecules were investigated in [H, H^-HJ , with the 
help of the coherent state method. The type of phase 
change discussed there, important for the context of stud- 
ies presented in this contribution, is related to the SU(3) 
and 50(4) limit and the transition between them in a 
molecule, which can consist of two atoms or of two nu- 
clei. Here wc will restrict ourselves to nuclear clusters 
only. Two groups of models will be discussed: the Phe- 
nomenological Algebraic Cluster Model (PACM) and the 
Semimicroscopic Algebraic Cluster Model (SACM). The 
Vibron Model belongs to the group PACM [25|,[26|. In the 
PACM the minimal number of relative oscillation quanta 
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is always zero. In contrast, in the SACM there is a mini- 
mal number of relative oscillation quanta, no, due to the 
Wildermuth condition [23], which is necessary in order 
to observe the Pauli exclusion principle. Though there is 
a lot of investigation in the PACM on phase transitions, 
applications to actual nuclei are very rare. Onl y in [28| . 
from the beginning of the vibron model, and in [23 ] have 
we found applications to real cluster systems. 

Some of the important questions we would like to dis- 
cuss in this and a forthcoming publication are: What 
is the difference between taking into account or not the 
Pauli exclusion principle? What are the orders of the 
phase transitions in the models discussed? How docs 
one define the thermodynamical limit? Normally, only 
second-order phase transitions appear between the SU (3) 
and SO (4) dynamical symmetries @. So, is it also pos- 
sible to obtain, under certain circumstances, a first order 
phase transition? 

This contribution restricts itself to the geometric map- 
ping of an algebraic Hamiltonian within the PACM and 
SACM. Already there some important differences arise. 
One main result will be that in order that the PACM 
reproduces the same results as the SACM, higher order 
interaction terms are necessary which simulate the effects 
of the Pauli exclusion principle. Differences and common 
features between the PACM and SACM will be discussed. 
The main reason for the differences is the large overlap of 
the clusters, making it necessary to antisymmetrize the 
many nucleon system. The PACM, which ignores the 
Pauli exclusion principle, will consequently fail in satis- 
fying basic conditions. For atomic molecules this prob- 
lem does not arise, because the two atomic nuclei are 
separated in space, and thus no exchange effects play a 
role. Some caution must be exercised using the compari- 
son; because the structures of the individual clusters are 
described within the SACM by the shell model, we are 
obliged to compare to a PACM that also uses the SU (3) 
model. In general, the IBA model has been used in the 
literature [3, [29| . Additionally, the parameter by which 
the number operator of 7r bosons is multiplied is fixed in 
the SACM because it describes the mean field. Within 
the PACM this parameter can be chosen arbitrarily. 

This all refers to the first question. The others will be 
addressed in the second part, dedicated to the study of 
phase transitions. 

The paper is structured as follows: In section Ql] the 
Semimicroscopic Algebraic Cluster Model is revisited, in- 
troducing some novel features, including the definition of 
the Phenomenological Algebraic Cluster Model (PACM) 
as a special limit of the SACM not observing the Pauli 
principle. In section UTT1 the coherent state formalism will 
be implemented in both models. The PACM coherent 
state can be recovered from the SACM when one sets 
the minimal number of relative oscillation quanta, no, to 
zero. In section ITVl the geometrically mapped potentials 
are derived for the two models, and finally, in section IVl 
conclusions are drawn and a discussion is presented on 
the differences of the PACM to the SACM. 



II. THE SEMIMICROSCOPIC ALGEBRAIC 
CLUSTER MODEL RECONSIDERED 



In this section we present a brief overview 
of the Semimicroscopic Algebraic Cluster Model 
(SACM) [H, [T| and introduce further amendments of 
it necessary for our study. These new elements appear 
in all three subsections. Previous applications of the 
SACM concerned describing the spectroscopic properties 
of core+a-type [30| and other [3l[ two-cluster systems. 



A. The group structure 

We start with reviewing the vibron model [25j , in which 
the relevant degrees of freedom are oscillations in the rel- 
ative motion of two structureless clusters in three dimen- 
sions. The operators describing them are boson creation 
and annihilation operators with angular momentum one: 



1 rn i 
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To this system one adds the spinless cr^ boson creation 
and rr annihilation operators. They define a cut-off, 
through the condition that the total number of bosons 
JV = n v + n<j is kept constant. The cr-bosons have no 
physical significance, which will play a role later on if 
one intends to define a thermodynamical limit. The rc m 
operators satisfy the relation 



TV — {-!) 7T_ m . 

The sixteen boson number conserving operators 
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act as the generators of the Ur(4) group, where R stands 
for relative motion. There are two subgroup chains that 
contain the SOr(3) rotation group. The irreducible rep- 
resentations of the subgroups supply quantum numbers 
to define bases that are associated with the two dynam- 
ical symmetries: 



U R (4) d 
[N, 0,0,0] 

where 
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n n = N, N — 1,...,1,0, 
Lr = n-n, n n — 2, 1 or 0, 



Mr = Lr, Lf> — 1 
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and 



U R (4) D SO r (4) D SO R (3) D SO r (2) 
[N, 0,0,0] (u,0) L R Mr, (6) 

where 



uj = N,N — 2, 1 or 0, 
Lr = uj,uj -1, 1,0, 

Mr = Lr, Lr — 1, —Lr. 



(7) 
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The SU(3) dynamical symmetry is generally believed to 
be the vibrational limit of the system around a spherical 
equilibrium shape, while the 50(4) dynamical symmetry 
describes static dipole deformation. 

The vibron model formalism reviewed up to this point 
handles only the relative motion of the clusters and ne- 
glects their internal sructure. In order to incorporate 
these degrees of freedom too, the SACM applies Elliott's 
SU(3) model [H|]. The orbital structure of the clusters is 
then described by the SUc k (3) group, where Ck refers to 
the k th cluster, k = 1, 2. The Elliott model applies LS 
coupling, but in many cases the 5 spin degree of freedom 
does not play a role. This is the case, for example with 
even-even clusters, and for the sake of simplicity we shall 
consider clusters of this type in what follows. 

It is essential that in the SACM the SU (3) group ap- 
pears not only in the description of the relative motion 
and the individual clusters, but also in the description of 
the unified nucleus. The typical group structure associ- 
ated with a two-cluster system in the SACM is then 

SUcM ® SU C2 (3) (g> SU R (3) D SU c {3) <8 SU R (3) D 

(Ai,/ii) (A 2 ,/xa) (n*,0) (Ac,Mc) 

SU(3) D 50(3) D 50(2) 

(A,//) kL M, (8) 

where (Afc, Hk) refers to the SUc k (3) irreducible represen- 
tations (irreps) of the individual clusters, which are then 
coupled to intermediate irrep (Ac,Mc)- These irreps are 
the ones associated with the ground-state configuration 
of the k th cluster, n, is the number of relative oscillator 
quanta, while (A, fi) is the total SU (3) irrep. L and M 
are the angular momentum and its projection, and k is 
used to distinguish multiple occurrences of a given L in 
(A, A*)- 

The model space of the SACM is obtained by compar- 
ing all possible irreps (A, /z), as given in ©, contained 
in the product (Ad,MCi) S3 (Ac 2 ,Atc 2 ) ® («7r>0), with 
those of the shell model and retaining only those irreps 
which appear in the shell model. Computer codes deter- 
mining the model space are available and can be obtained 
on request. In most cases, however, it is easy to retrieve 
the irreps by hand. In this manner the Pauli exclusion 
princi ple is observed (for some illustrative examples, see 
Refs. [H Gj|). The SU{3) basis is also useful in elimi- 
nating the spurious center of mass motion. 

We note that the above 5i7(3) matching procedure also 
reproduces the Wildcrmuth condition [27 1 in a natural 
way. This condition prescribes a minimal number of os- 
cillator quanta (i.e. n w ) in the relative motion. Apart 
from the case of closed-shell clusters, however, it is only 
a necessary condition for the handling of the Pauli exclu- 
sion principle. 

It is now worthwhile to discuss the possible dynamical 
symmetries of the SACM based on those of the vibron 
model. The SU (3) dynamical symmetry is clearly asso- 
ciated with the (O group chain. The equivalent of the 
50(4) dynamical symmetry of the vibron model, how- 



ever, can be considered only an approximate dynamical 
symmetry in the SACM. The reason is that due to the 
Pauli principle part of the set of 50(4) basis states has to 
be excluded from the model space. Although is not a 
good quantum number in the 50(4) limit, the 50(4) ba- 
sis states can be written as linear combinations of SU (3) 
states, so excluding these below the minimal allowed 
value distorts the structure of the 50(4) basis. Finally, a 
third dynamical symmetry can also be derived from the 
SU(3) dynamical symmetry of the vibron model. The 
group structure associated with this 50(3) dynamical 
symmetry is 



SU C (3) ® U R (4) d 50 c (3) ® 50^(3) D 50(3) D SO (2) 
(X c .fxc)[N, 0,0,0] L c L R L M. 

(9) 



The difference between the previously mostly ignored 
© chain and the one appearing in (J5J is of dynamical 
nature in the sense that the interaction in the former 
case does not contain terms typical of the coupled SU(3) 
degrees of freedom. In fact, the SU(3) groups do not play 
a role other than supplying labels for classification of the 
states. In terms of interactions we can call the scenarios 
associated with the © and ([5]) chains as weak and strong 
coupling limits, respectively. The two limits are the same 
when the two clusters are both closed-shell nuclei, but 
when at least one of them is not (i.e. its internal (A&, /ik) 
irrep is different from (0,0)), a clear difference between 
the two limits arises. 

Before closing this subsection it is worthwhile to com- 
ment on the typical selection rules characterizing the dy- 
namical symmetries. This is also related to the band 
structure determined by the appropriate group structure. 
In the basis associated with the SU(3) dynamical sym- 
metry of the SACM the bands are defined by the (A, /i) 
and k quantum numbers (see Eq. ([5])). where k is ob- 
solete when either A or fi is zero, as is the case in the 
SU(3) limit of the vibron model too (see Eq. The 
states belonging to the same SU(3) irrep are connected 
by the quadrupole operator, the SU(3) tensorial charac- 
ter of which is (1,1). This operator leaves and the 
parity intact and changes the angular momentum by two 
units, so it describes electric quadrupole transitions. On 
the other hand, bands associated with the 50(4) dy- 
namical symmetry are characterized by the w quantum 
number defining the 50(4) irreps (sec Eq. ([5])) and con- 
tain states with both even and odd angular momentum, 
i.e. with both positive and negative parity. The in-band 
transitions are described by the 50(4) generators, which 
play the role of the electric dipole operator. The two dy- 
namical symmetries thus lead to different selection rules, 
and this has to be taken into account when they are ap- 
plied to some concrete physical problem. 
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B. The Hamiltonian 

Let us now turn to the Hamitonian associated with the 
SACM. While in most typical applications it is sufficient 
to consider interaction terms constructed as two-body 
terms, here we argue that a specific third-order inter- 
action term is also necessary to stabilize the spectrum. 
Furthermore, as another new element we shall separate 
the Hamiltonian into terms associated with the three dy- 
namical symmetries identified above. The parametriza- 
tion introduced this way allows interpolation between 
the dynamical symmetries, changing certain parameters 
of the Hamiltonian, like x and y (see equation below). 
We consider two cases: i) both clusters spherical and 
ii) one spherical cluster plus a deformed one. Exam- 
ples for these two scenarios are the 16 0+a — > 20 Ne and 
20 Ne+a — > 24 Mg systems, examined in Paper-II of this 
series. 

The Hamiltonian is given by 
H = xyH su(3} + 2/(1 - x)H so{i) + (1 - y)H 0(3 ) (10) 
with 

H S u (3) = hurt* + a C ius C 2 (A c , "c) 

+(a- bAn- K )C 2 (A, fi) + (a - bAn^)C 2 (n OT , 0) 
+ 1 L 2 + tK 2 

it T 2 i (1) T 2 

-"SO(4) = 0, C -Lc +0-R 

+ 7j L 2 + ^[( 7 rt. 7r t ) -( a t ) 2] [ (ir . ff) _( ff) a] 

H (3) = hujn^ 

+a c L c 2 + a^L 2 R + 1 L 2 , (11) 

with An^ = tin — no, no being the minimal number of 
quanta. The acius is the strength of the quadrupolc- 
quadrupole interaction, restricted to the cluster part, 
while R and C denote the contributions related to the 
relative and cluster part respectively. Further interac- 
tion terms are the total angular momentum operator, L 2 , 
and the K 2 operator, defined in [TH, [l9| which classifies 
the rotational bands, giving the projection of the angu- 
lar momentum onto the intrinsic z axis. For the case of 
two spherical clusters, the second-order Casimir opera- 
tor of SU{3) is just given by ^(n^ + 3). Note that in 
the case of deformed clusters the information about the 
deformation only enters in the SU(3) dynamical limit. 

Note that the division in (jlip is done according to dy- 
namical symmetry limits and not according to two terms, 
referring to each cluster, one to the relative motion and 
one to the interactions between them. If one wishes to do 
that, all what has to be done is to decouple the different 
contributions. For example, the L 2 operator can be writ- 
ten as (L c + L R ) 2 = [L 2 c + L 2 R + 2 (L c ■ L R )] . The 
first and second term refer to the cluster and relative an- 
gular momentum, respectively, while the last term refers 
to the coupling between the channels. This can be fur- 
ther divided by writing the cluster angular momentum as 

Lc-m = {Ldm + Lc 2 m)- 



The division according to dynamical symmetries in 
(fTTj) was done as follows: In the SU(3) limit the cou- 
pling of interaction operators is on the level of £[7(3), 
i.e., Casimir operators of SUc(3) and SU(3) have to ap- 
pear. This is called the strong coupling limit. The 5*0(3) 
limit couples only at the 50(3) level, i.e. only the angu- 
lar momentum operators appear (no second- and higher- 
order SU(3) Casimir operators, except n T and functions 
in it). This is called the weak coupling limit. This is 
reflected by the appearance of interaction terms related 
only to angular momentum (L R , L 2 -, and L 2 ). In the 
literature one usually refers to this latter limit as the 
SU (3) limit, where the Hamiltonian contains only the 
fkurin term plus at most some weak anharmonic terms. 
Here we feel it necessary to change the notation because 
we understand the SU(3) limit to include also terms such 
as a strong quadrupole-quadrupolc interaction. In order 
to compare with results in the literature, this has to be 
kept in mind when we report on phase transitions be- 
tween different dynamical symmetries. The 5*0(4) limit 
is defined through the appearance of the second-order 
Casimir operator of 50(4). This limit is called the de- 
formed limit because the interaction will always produce 
a potential with a deformed minimum. 

In principle, one can add the angular momentum oper- 
ator of the deformed clusters (L 2 , k = 1,2). We exclude 
this interaction for the moment. 

The new higher-order interaction appearing in the 
third term of H su(3) needs some explanation. The whole 
term is related to the quadrupole-quadrupolc interaction, 
which is present in any nuclear system. However, with- 
out the —bAn-K (with —6 > 0) correction, states which 
contain a sufficiently large n w will be lower in energy 
than states with the minimal number of tv bosons, no. 
This is due to the dependence on n 2 in the second-order 
Casimir operator, which will finally dominate over the 
Tiuin^ term for a sufficiently large number of 7r-bosons. 
In the standard treatment, when n^ is conserved, a sim- 
ple restriction to small An w suffices to circumvent the 
problem, i.e., states with large An^ are simply not taken 
into account in the model space. 

This effect was studied in [33| within the context of 
the symplectic model of the nucleus [l^, H3, IH|. Also 
there, the quadrupole-quadrupolc interaction dominates 
over the kinetic energy and finally will promote high n w 
states to low energies, even below the physical ground 
state. This problem was solved by subtracting from the 
quadrupole-quadrupole interaction the so called Trace 
Equivalent part [33[ , which insures that the average mean 
field is still represented by a harmonic potential. When 
no correction is applied, the mean field shell structure 
is destroyed and a mean oscillator structure, one of the 
main assumptions of the shell model, cannot be assumed 
anymore. This was also noted within the SACM in fl9| . 
where correction terms of the type An T , mentioned here, 
were included. Without these corrections the problem in- 
creases significantly when interactions mixing states with 
different n T arc considered. Then, avoiding states with 
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large n„ is not an option, as is in the case of conserved 
n w , when the model space can be limited in using 
physical arguments. 



PACM: the phenomenological limit of the 
SACM 



The minimal number of w bosons is an essential re- 
quirement in the SACM to incorporate the Pauli princi- 
ple. However, the formalism allows setting this minimal 
number to zero. This limit of the SACM can be de- 
fined as the Phenomenological Algebraic Cluster Model 
(PACM). It has to be stressed that the difference be- 
tween the SACM and the PACM manifests itself only in 
the model space, while the two models share the same 
Hamiltonian and other operators. Obviously, the differ- 
ent model space will lead to different matrix elements in 
the two models. Note that the minimal number of rela- 
tive oscillation quanta is either or no. It is not allowed 
to choose a number in between, because each of such a 
number violates the Pauli exclusion principle. 

When both clusters are closed-shell nuclei, the PACM 
essentially recovers the vibron model (25[. 

One of the main objectives of the present work is to 
investigate the similarities and differences between the 
two approaches. This is especially interesting within the 
context of the coherents state formalism, because in other 
models restrictions similar to those in the SACM (i.e. re- 
stricting the boson number) are unknown. In this sense 
the formalism of the PACM is closer to that of other 
models. Due to the minimal number of ir bosons the 
formalism of the SACM will obviously become more in- 
volved. It is our aim to explore this conflict between 
the physical importance of a fundamental principle (i.e. 
the Pauli principle) and the technically more complicated 
formalism arising due to it. 



III. COHERENT STATES AND THE 
GEOMETRICAL MAPPING 



In this section the coherent state is presented, which is 
used to obtain a geometrical mapping of the SACM and 
PACM in the next section. 

The use of coherent states is the most common method 
of applying a geometrical mapping 0, d, 0, @, HO, [13 • One 
advantage is that the coherent state can be expanded 
in terms of the complete set of states for a given total 
number of bosons, N, (in the SACM, this refers to all 
allowed basis states for a given total number of quanta). 
The ground-state energy is usually reproduced very well. 
The coherent state also provides a transparent relation 
to collective variables. Its use is justified by noting that 
it corresponds to the Gaussian Overlap Approximation 
within the Generator Coordinate Method [36| . skip ping 
the term of the zero-point motion. As shown in |3g|, this 
method allows the definition of a potential with usually 



good results. However, the mass parameters of the ki- 
netic energy arc usually not reproduced very well. In 
order to obtain a kinetic energy too, the coherent state 
variables have to be defined as complex variables 0, 0] ■ 
We do not consider the kinetic energy due to the reason 
mentioned above, and focus on the potential. 

The coherent state within the SACM was introduced 



17] 



M N , no {a-Tt) na [ < r t + (a-7rt)] 7V |0) 

N\ d™° r + , ^N+n 



N,n 



(7V + n )!d 7l io 



7i (a ■ tt^)]' 



10) , 
(12) 



where, for convenience, we redefined the total number of 
relative oscillation quanta as (N + no), while the 71 pa- 
rameter has to be set equal to 1 after the differentiation. 
The normalization factor is given by (T7| 



iV! 2 



N, 



(yV + no)!d 7 r° d 7: 



^ [1 + 7i72(a* • a 



iiV+no 



(13) 

Again, the 7^ have to be set equal to 1 after the applica- 
tion of the derivatives. 

The a is a short-hand notation for the, in general, 
complex variables a m (m = 1, 0, —1). The coherent state 
with complex a coefficients is the most general linear 
combination of the boson creation operators. For static 
problems the requirement @, 0] 



(14) 



reduces the number of real parameters to three, namely 
to ao and the real plus the imaginary part of aq_i. 

In the Appendix we present the results for the ge- 
ometrical mapping for the important interaction terms 
appearing in the Hamiltonians. We define 



(a 



a 



-1 



a m a- 



(15) 



where a represents a measure of the inter-clustcr distance 
[ItJ and a 2 a short-hand notation for (a ■ a). Because 
the only relevant variable is the inter-cluster distance, we 
can express the potential in terms of this sole variable a. 

The importance of the coherent states resides in the 
fact that they provide us with the possibility to define a 
Potential Energy Surface (PES) 

V{a) = (a\H\a) , (16) 

in terms of the collective variables a m . 

A. Renormalization of the variable a 

It is often convenient to transform a to other re- 
lated variables. In the literature there are differ- 
ent conventions, which often contradict each other 
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|22h24 l37H39j . Here we present some arguments to 
justify our choice, restricting for simplicity to the PACM, 
where no = 0. 

In 0, [?J the renormalization of the interaction param- 
eters is proposed using the following reasoning: the ex- 
pectation value of a one-body interaction with respect to 
the coherent state is proportional to the total number of 
quanta. As an example, the expectation value of ain^ is 
given by 



(ain n ) = a±N 



(1 + a 2 ) 



(17) 



Assuming that a is of the order of one, this expectation 
value increases with N, which is unnatural because (n w ) 
is of the order of one (note that this argument assumes 
that a is of the order of one, too). In order to avoid this, 
it is recommended in [f| @] to redefine the interaction as 



N 



(18) 



In this manner, it is expected that the parameter a[ does 
not change significantly with increasing N. A similar ar- 
gument holds for the two-body interaction, recommend- 
ing to divide the corresponding interaction parameter by 
N(N - 1), and so on. 

This argument is supported by considerations given in 
[icj . There it is shown that one- and two-body inter- 
actions scale like N and N(N — 1), respectively. This 
dependence has to be canceled in order to define a ther- 
modynamic limit N — > oo. However, the N used in (40| 
can always be related to the number of particles. For 
example, in the IBA [U N is given by half the number 
of valence nucleons and in Lipkin-type models (two-level 
systems) the N is given by the number of states in the 
lower level which are completely filled in its lowest states. 
In the algebraic cluster models, however, the iV is not re- 
lated to a number of particles but rather to a boson cut- 
off. The er-boson is introduced merely to define a cut-off 
for the number of 7r-bosons (n„ < N). Physics requires 
that the final results do not depend on this cut-off. 

In what follows we present a proposal on how to treat 
the cut-off in algebraic cluster models. Note that in (TT5)) 
a conceptual problem arises already, which has not been 
considered before in the literature. Here we try to con- 
vince the reader that the former renormalization (simply 
dividing the interaction parameters by a power in N) is 
too simple and a more sophisticated renormalization has 
to be applied. Here we illustrate the situation in three 
different ways. 

a) Consider first the hourly term, which represents the 
mean field of the harmonic oscillator in the SACM. Di- 
viding it by N implies that Two 1 = ^ tends to zero as 
N — > oo. The contradiction becomes apparent because 
fiuj is a physical value, which has to be kept independent 
of N, while a\ in Eq. (|17[) is used as an abstract param- 
eter. 

b) Looking at it from a different angle, let us consider the 



operator n n alone. When applying in an harmonic 
oscillator basis, in the lowest states the eigenvalues of 
n T are small numbers. The geometrical mapping gives 

2 

N ^° a ^ , which suggests an increase proportional to N. 
The only way to maintain numbers of the order of one is 
to redefine a by 



5 1 
N 



(19) 



such that the a tends to zero for N — > oo. In this way, 
with 8 of the order of one, the geometrical mapping also 
gives results of the order of one. Thus, restricting for 
example to the PACM (n — 0), the coherent state in 
(fT2"]) is redefined as 



\S) = 



1 



N 



The expectation value of n T then leads to 



«77 = 



1 



N ) 



|0) . 

(20) 

(21) 

oo, 



The global dependence on N vanishes, while for N 
the TV-dependence disappears. 

The same happens for most of the two-body interac- 
tions, as we will discuss further below, with some differ- 
ences when the er operators are involved, 
c) The above choice of the new variable S can be justified 
in a third way. To prove this, let us consider the coor- 
dinate operator defined without er-bosons, i.e. without 
cut-off (N), namely 



h 



2mu 



(22) 



When the cr-bosons are introduced, this operator has to 
be changed. The new operator, called the algebraic co- 
ordinate operator, should satisfy the following minimal 
conditions: i) The total number of bosons has to be kept 
constant, i.e., each tz^ has to be multiplied by er and 
each 7r m has to be multiplied by a cr; ii) the definition 
of the distance operator should be independent of the ba- 
sis and Hamiltonian used; and hi) for N — > oo it should 
converge to the standard form given in (|22|) . The pro- 
posed algebraic coordinate operator is given by 



2Nmuj 



(T^TY, 



(23) 



The "a" refers to an algebraic operator. The operator 
itself does not change the total number of bosons, as 
required by the above condition i). The N in the de- 
nominator of the square root is introduced because in 
the harmonic oscillator basis the matrix elements of the 
CT-operators behave like y/N — n w , which for large N and 
small number of tt bosons is approximated by y/~N. This 
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approximate value of the cr operators is satisfied in any 
basis, with the condition that the average number of ir 
bosons is much smaller than N (though, the structure 
is particularly simple in the harmonic oscillator basis). 
Thus the 1 / y/~N factor cancels approximately the contri- 
butions due to the addition of the er^ and cr operators. 
In this form, the algebraic coordinate operator does not 
depend on the basis used (the can be applied to any 
kind of basis) and nor on the Hamiltonian, thus, satisfy- 
ing condition ii) . For very large N the expressions of the 
physical and the algebraic coordinate operators tend to 
each other, satisfying condition hi). 

This definition agrees with (37l . l38j where an algebraic 
model for atomic molecules is discussed. Often (see for 
example [39[) one defines the radial distance as a func- 
tion of the dynamical symmetry, relating it indirectly to 
the matrix element of the dipolc operator, without any 
further considerations. This violates condition ii) above. 
We insist that the definition of the radial coordinate op- 
erator has to be independent of the Hamiltonian in the 
Hilbert space considered. The Hamiltonian determines if 
there is a dynamical symmetry or not, which should be 
independent of the radial coordinate operator, while the 
Hilbert space as such is independent of the basis used. 

The expectation value of the algebraic coordinate op- 
erator is 



2NH 



mix) (1 + a 2 ) 



(24) 



We define this as the algebraic distance r" , which is, by 



definition, of the order of one. 
gives 



Inverting this relation 



(l + a 2 ) 



2NH 



(25) 



which again provides the dependence of a m on N. It 
suggests we redefine a m in terms of S m and N as given 
above, i.e., 



giving 



{( 



- r 



a\2 
0/ 



N > 



N 



2muj 



2muj J i _l §_ 

N 

(28) 



where the arrow gives the limit for large N. As long as 
the expectation value of the algebraic distance operator 
is greater than the square root of this expression, it is 
safe to identify the rg as the distance between the two 
clusters. The square root of (|28[l gives numbers of the 
order of 1 fin. 

When applying the PACM the usage of another vari- 
able 



(29) 



(l + a 2 ) ■ 



will be found more favorable. The range of 1 is 
< /3 2 < 1, because the range of a is < a 2 < oo. 
The /3 2 variable can be related to S 2 as it was to a 2 . The 
main reason to use /? is that within the PACM it simpli- 
fies most expressions, as can be seen further below. For 
the SACM we will return, for convenience, to the variable 
a. 

Before closing this subsection we note that in the 
SACM the relation between a and the relative distance 
of the two clusters is given by a ~ (r — ro), fLl\ , with 
r o ~ y/no- This means that there is a minimal differ- 
ence in the inter-cluster distance due to the presence of a 
minimal number of 7r-bosons, i.e., clusters can not over- 
lap completely. In the PACM, no = and thus ro = 0. 
Defining the minimum at a = as the "spherical" min- 
imum then loses its meaning, because in the SACM it 
already corresponds to a minimal distance. Nevertheless, 
we will continue to call a minimum at a = a "spherical" 
minimum and a ^ as a "deformed" minimum. 



JV 



(26) 



IV. THE GEOMETRICALLY MAPPED 
POTENTIAL 



with S m given by 



2h 



(27) 



a dimcnsionless measure of the distance between the two 
nuclear clusters. We claim that this is a consistent way to 
define the distance between clusters trough the variable 



As discussed previously, the difference between the 
SACM and PACM model spaces manifests itself in the 
difference of the matrix elements of physical operators, 
even if the operators themselves are the same in the two 
approaches. In this section we determine the potential 
energy surfaces in both models and explore the relation 
between them. 



The validity of (|24[) . the definition of r"„, depends on 
the fluctuations of the related expectation value. This 
will be discussed shortly. The inter-cluster distance vec- 
tor can always be chosen along the z-axis. The square of 



the variation ((rg 



<(2\ 



can then also be calculated, 



The SACM case: Pauli principle taken into 
account 



Applying the coherent state for the SACM to the 
Hamiltonian (|I0[) , one obtains the geometrically mapped 
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potential 



(H) 



C(x, y)-(b+ b)xy A(x, y)a 



2 F 11 (a 2 ) 
F 00 (« 2 ) 



-B(x,y)c 



F 22 (a 2 



6^33 K 



F 0Q (a 2 ) F 00 (a 2 ) 



-C(x,y)a 



2 x 20 



^oo (« 2 ) 



(30) 



where 



C(x,y) 



((actus + a + bn ) xyC 2 (Ac, Mc) + + (1 - xy) L c 2 ) 

+ xyt(K 2 ) + ^(N + n ) (N + n -l)y(\-x) (31) 



and the Fij(a 2 ) functions are defined as 



F m a') = 



N+n 

x E 

fc-TlQ 



Fn (a z ) = 



(TV!) 2 



{N + n )l 





' fc! ' 


( w r°) 


[(k-n )\\ 


(Nlf 




(N + n - 


-1)! 



, 2k 



N+n -l 

x E 

fc— max (no — 1,0) 



+ n - 1 
fc 



(k + l)\ 



-i 2 



(fc + l-no)! 



F 22 (a 2 



(N\) 



(N + n -2)\ 

N+n -2 

x E 

k— max (no —2,0) 



N + n - 2 
fc 



(fc + 2) 



(fc + 2 - n ) 



Q 



2/.- 



<- 2 k 



JV+n c ,-2 

x E 



fc— no 





(TV!) 2 




= 


(TV + n - 2) 




/ N + n - 2 \ 


' k\ ' 


V 


* ; 


.(A -no)! J 


X 


' (fc + 2)! " 


a 2k 




_(k + 2-n Q )\_ 



F 33 (a 2 



(N 



(N + n -3)! 

x E ( ? 

k— max (no —3,0) 



(fc + 3)! 



(fc + 3-n )! 



.2fc 



a 



(32) 



Further, the constants appearing in (|30[) are defined as 
1 



A(x,y) = - 



(b + b) xy 



(fiu [yx + 1 — y] 



7 + (1 — yx)a 



B(x,y) 



C(x,y) 



+xy [a - b] [4 + r x + T 2 ] 

+4xy [a — b ] + xybno [4 + Ti + r 2 ] + Axynob 
-bxyC 2 (A c , u c ) - |y (1 - a:) (N + n - 1)) 

1 xy [a + a — 66 — 66 



(33) 



(b + b) xy 

>{ri + i 

fy(l-x) 



-fo{ri + r 2 } + n (6 + 6)] + -y(l-z) 



(b + b) xy 
where r^, according to fl7j . is given by 

r fc = ((Afc.^Jlo^-^^KAfc,^)) 



-TVo,fe/3 fc 

7T 



This was obtained b> 



a 2m (k) 

(34) 

eometric mapping of the sym- 



med by as 

plectic model [Io|, [3 El). The (Afc,/ifc) denotes the 
SU (3) irrep of the deformed cluster number fc. In the 
case that it is spherical, = 0. The N k is the sum 
of the total number of quanta n^ of the deformed cluster 
plus \(Au — 1), where Ak is the number of nucleons in the 
fc cluster. This last term is the zero-point energy with 
the contribution of the center of mass already extracted. 
The a 2m (k) is the deformation variable of cluster number 
fc. In (1341) we used only the m = component of a 2m (k) 
and defined it /3k, the deformation of cluster fc (not to 
be confused with the (3 variable appearing in Eq. (|2"9|l ). 
This implies that the deformed cluster is assumed to be 
axially symmetric and it is in line with the inter-cluster 
z-axis, which connects both clusters. When the z-axis 
of the deformed cluster is inclined with respect to the 
molecular z-axis by an angle 6, the deformation value 
fik is multiplied by a matrix element of the rotation ma- 
trix, which only changes the numerical value of ftk, i.e., 
Pk = ^oo(^)/^fe- For simplicity we do not include these 
orientations in the discussion. Furthermore, it will not 
change the basic results. 

In discussing the phase transitions it is possible to 
choose as the independent parameters of the theory A, 
B and C, which themselves are the functions of all in- 
teraction parameters of the theory. This structure will 
be used in the second paper, investigating the possible 
phase transitions and the phase diagram. 

For x = in ([10]), i.e. in the case of the 50(4) to SO(3) 
phase transition, the discussion has to be modified due 
to the xy factor appearing in the denominators in (|33l) . 
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In this case the potential maps to 



The relevant formulas arc 



V = (H) 
Aa 2 



n F u (a 2 



+C 



F 00 (a 2 ) 



Ba 



4 F 22 {a 2 ) 
F 00 (a 2 ) 



Ca 



rpN-2 



^oo (a 2 ) 



(35) 



with 



A = [hu[l-y] + 2 



7 + a« 



-v(N + n -l) 



B 

C 



2 ! > 



-B 



(36) 
(37) 

(38) 



In this case only two independent parameters A and B 
appear. 

Note that for c > 0, the C is positive (remember that 
— b > 0). When C < 0, the situation corresponds in the 
50(4) limit to a ground state where all bosons are de- 
coupled and the highest state is the one where all bosons 
are coupled in pairs. 

In (f3"0"|) we have to add a constant term, such that the 
geometrically mapped potential is zero at a = 0. This 
is a permitted renormalization of the zero-point energy. 
This constant will be determined further below. 

It also has to be noted that the 50(4) dynamical sym- 
metry needs special care due to the truncation of the 
SU (3) basis required by the Pauli principle. Eliminating 
these components leads to the destruction of the 50(4) 
dynamical symmetry. However, we will still denote it a 
50(4) dynamical symmetry, because the operators in the 
Hamiltonian will be the same. Instead of the 50(4) basis 
we will stay within the SU(3) basis, because only there 
can the Pauli principle be implemented easily, canceling 
all states with < uq. Note also that the total number 
of bosons is now (N + no) and not just N. 

A very useful consideration is the investigation of the 
potential in the a — > oo and a — > limits. In the first 
limit we will see that the potential approaches a constant 
value depending on (N + no), which is due to the finite 
size of the boson space. For large values of a the coherent 
state contains only 7r-bosons and cannot increase the en- 
ergy any further. The second limit (a — > 0) is necessary 
to adjust V(a = 0) = 0. These consideration will be im- 
portant in the second paper, when the general structure 
of the SACM phase diagram will be discussed. 

i) Limit a —> oo: 



2^11 

Y 

Fqo 

4 F22 

Y — — 

Fqo 

,6^33 

Fqo 

N-2 



TP 
2 r 20 



(N + n ) 

(TV + noMiV + no-1) 

(N + n ) (N + n - 1) (N + n - 2) 
1 



F 



0(1 



N(N-1) 







(39) 



With this, the limit of the complete geometric potential 
(I3T1) is given by 

V -> C - (b + b) xy {A (N + n ) - B {N + n ) (TV + n - 1) 
+ {N + n ){N + n -l)(N + n Q -2)} . (40) 

Depending on the signs and values of A and B, this limit 
is either positive or negative. For the positive value the 
limit for (N + hq) — > 00 is then +00, leading to a stable 
potential, while if it is negative the limit leads to —00, 
leading to an unstable potential. 

ii) Limit a — > 0: 

The relevant formulas are 



2-^11 

^00 
4^22 

6 F 33 

a ~F~ 

-TOO 

F N-2 
,2 ^20 



Ft), 



n 




n (n 


-1) 


n (n 


- 1) (no 


7V(7V- 


- 1) (no 



(41) 



With this, the limit of the complete geometric potential 
(|5T1) is given by 

V(a = 0) -> C - (b + b) xy {An - Bn (n - 1) 

+n (n - 1) (n - 2)} (42) 

which is independent of N. This result can be used to 
adjust the potential to zero at a = 0. 



B. The PACM case: Pauli principle not taken into 
account 



Using the Hamiltonian as introduced in the section HTl 
and the coherent state of section IIII1 for the case when 
the Pauli exclusion principle is not taken into account, 
the potential is obtained by calculating the expectation 
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value of the Hamiltonian as 



(H) = Vtf) 



(actus + a) xyC 2 (A c , fi c ) + -N(N-l)y(l- x) 



Nf3 x 



(Twj (xy+l-y) + (a- b) (4 + Ti + T 2 ) 



-4 a 



-bC 2 (A c ,Mc)) -(l-x)yc(N -1) 
+2 (j+(l-xy)a^ 



+ N (N - 1) (3 4 [xy (a + a - 66 - 6b - b (Ti + T 2 )) 
+ (1 - x)yc] 

- N (N - 1) (N - 2) /3 6 a;y (6 + fe) + C 2 (A c , He) axy 



+ - (l-x)ycN(N-l) 



(43) 



Denning 



.4 



B 



-[(b + b)xy(N-l)(N-2)] 1 
x \hcj (yx + 1 - y) + 2 (7 + (1 - 2/x)a^ 1} ) 
+ 4.xy (a - 6) + xy (a - 6) (4 + Ti + T 2 ) 
- bxyC 2 (A c , Mc) - 2/ (1 - x) c (N - 1) 
xy(a + a-6(b + b) -b(T 1 +T 2 )) +cy(l-x) 



(N —2) (b + b) xy 
C = ({acius + a)xyC 2 (\ c ,Hc) +-fL c 2 + (1 - xy) L c 2 ) 

(44) 



+xyt(K 2 ) + C -N(N -l)y{l-x) 



the potential acquires the form [4lj, |42J 

V = N(N-l)(N-2)(-(b + b)xy){Af3 2 -Bf3 4 + f3 e } 
+C , (45) 

which allows us to define a new, normalized potential 

V = {A(3 2 ~ Bf3 4 + f3 6 } . (46) 

In the definition of V wc extracted the factor (b + b)xy, 
such that there appears no factor in front of the ft 6 term. 
This poses no problem as long as x is varied from to 1 . 
In the limit of x the A and B values also approach 
±00, depending on the sign. However, for the 50(4) to 
50(3) transition, the x value is always zero. For this case 
we include the x value within the parenthesis, yielding a 
vanishing factor of the sextic term. 

Comparing the potentials obtained from the same 
Hamiltonian in the SACM and PACM approaches leads 
to a remarkable finding. The potential in the SACM 
framework is rather different from its PACM counterpart, 
however, a similar potential can also be generated within 
the latter framework too. This can be achieved by includ- 
ing higher-order interactions of the type F^n^) / F^n^) , 
with appropriate functions Fkin^). This demonstrates 
that observing the Pauli exclusion principle acts as if one 



used high-order interactions in a model which does not 
observe the Pauli exclusion principle. In fact, the non- 
linear terms simulate the presence of the Pauli exclusion 
principle. 



V. CONCLUSIONS 

In order to investigate possible phase transitions be- 
tween different limits corresponding to various dynamical 
symmetries, we reparametrized the Hamiltonian of the 
Semimicroscopic Algebraic Cluster Model (SACM) such 
that it allowed interpolation between the three possible 
limits. These were the strong coupling limit (SU(3)), 
the deformed limit (SO (A)) and the weak coupling limit 
(SO(3)). The latter limit was proposed in the present 
work and it differs from the strong coupling limit in the 
level on which the interaction terms of the relative motion 
and those of the internal cluster structure are coupled: in 
the weak coupling limit this is done on the SO(3) (i.e. 
angular momentum) level, while in the strong coupling 
limit the SU (3) algebra plays a role, introducing, e.g. 
quadrupole-quadrupole interaction between the two sec- 
tors. In the case of a system with two spherical clusters 
the weak coupling Hamiltonian is a simplified version of 
the strong coupling one, so it does not stand as a separate 
limit in itself. The SO(A) limit also has its limitations 
due to the truncation of the model space in the quan- 
tum number. 

The Phenomenological Algebraic Cluster Model 
(PACM) was introduced as a special limit of the SACM 
with the minimal number of the 7r bosons set to zero. 
This choice corresponds to neglecting the effects of the 
Pauli exclusion principle. Although this means giving up 
a fundamental physical requirement, this decision was 
inspired by the fact that the formalism of the PACM is 
closer to other similar models using the coherent state 
method. It appears instructive to study the differences 
and similarities between the SACM and PACM within 
this latter approach. 

The present work is meant to be the basis for a further 
study in which phase transitions are investigated by inter- 
polating between two dynamical symmetry limits. This 
method requires the application of large boson numbers, 
so as another new ingredient, the Hamiltonian was imple- 
mented with a third-order term in order to stabilize the 
energy spectrum in this situation. The potential energy 
surface was constructed in terms of a variable controlling 
the relative distance of the clusters. This was done both 
in the SACM and PACM framework. It was found that 
the potential obtained from the SACM can be reproduced 
within the PACM approach too by including higher-order 
terms in the Hamiltonian. This indicates that studying 
only the Hamiltonian, the effects of the Pauli principle 
can be simulated by higher-order interactions. 

The present results will be used in a forthcoming pub- 
lication that focuses on phase transitions beween phases 
determined by different dynamical symmetries of the 
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SACM and the PACM. 



is a real number. 



APPENDIX A: The Coherent State for the SACM 

We choose the most general structure for the coher- 
ent state, allowing arbitrary parameters, a m , which only 
coincide with a m when the static problem is considered. 
This will be important in future work, when we intend 
to treat the cranking formalism within the PACM and 
SACM, similar to the formalism presented in [43|, |44j . 
Nevertheless, as long as we are only interested in the po- 
tential energy surface for systems without rotation, the 
parameters a m will form a simple tensor. 

We use the definition 



(a. ■ tt') = 



7 / a m^ m 



(47) 



The formulas are now similar to those of Ref . [13] , but 
without the use of a possible tensor character of the a m 
and with the appearance of complex conjugate a* m and 
a m . One of the interesting matrix element is given by 

[3 



-to i [ s i 



) = (N + n Q )[a* xq]^ 
d"° d"° 

A no A no 7172 [1 + 7l72(<** ' ") 

d 7l d7 2 



iV! 2 
10 (N + n )\ 



N+n -l 



The a m are in general complex and arbitrary. The com- 
plex conjugate is denoted by a^. We also use 



(48) 



This will be important when we apply 



right. 



to the coherent state on the 



The conjugate coherent state is given by 



(a\ = Mi 



[er + (a* ■ 7r)] (a* ■ tv 



Nn ( 



with 



(a* -it) = 



E a *n* T 



, (49) 



(50) 



Thus, the ir^ acts on the left as an annihilation operator. 
Note that here we do not assume a tensorial behavior of 
the a m , contary to what we used in the body of the 
paper. In order to relate this a m to the one used in the 
paper, we have to assume a* m = (—1) ct- m - This is 
justified for a static problem, as discussed in the paper. 
The situation changes, when for example the cranking 
formalism is applied or not only the potential is intcnted 
to derive but also the kinetic energy. 

Mnti is the normalization factor, given by 



Nl 2 



AJlo AJlo 



(N ■ 



n )!d 7 r 



d72 



[1 + 7172(0;* • a 



(51) 

1, after performing the derivation, 
to the left, commutators of the type 
[•7r mi , 7r„ l2 ] appear and will give expressions proportional 
to a* m . However, acting with iv m , as it appears in a cou- 
pled expression, to the right, it will give expressions pro- 



taken at 7! = 72 = 
Acting with n\ 



portional to a m = (— 1) 



because we have first 



to lift the index of the annihilation operator, obtaining 
(-l) 1 -™-*-" 1 . Note that 



(a* 



a 



E 



(52) 



We have 



[«* x aim 



(lmi, lm2\Sm)a^ ll a ri 



(53) 



where the coupling sign "x" instead of "(g)" was used in 
order to indicate that we do not couple tensors. This is 
just a short-hand notation. 

In the geometrical mapping one has to take into ac- 
count that 



(a ■ 7r T ) 



(54) 



and thus 



[7r m , (a ■ 7^)] = (-1) 1 m a- m = a m 

[(a* ■ n),Tv} n ] = a* m . (55) 



Note the difference in the phase. 

In this sense, the mapping of different operators is com- 
pletely parallel to the one given in [I?} , with the exception 
of the definition in the coupling of a m . The geometrical 
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mapping of more relevant operators is given by 
» m | *-* / 

(Nl) 2 



a\cr'TT m \a) = 

-l) 1 - m (N + n )a m Af 2 Nno 

AJi AJio 



(iV + no)! 
72[l+7i72(a*-a)] Ar+ "°- 1 



1—771 



* ^ (W) 2 



(N + n )\ 



d"° d"° 



72 [1 + 7i72 (a* ■ a)] 



N+n -l 



a\a^cr\a) = N 2 



A/" 2 

k^Trtl [Sl] ® \n ® ir] lS2] 



[S3 



I") 



iV + 7i )(iV + 7i -l)A^ no 



2 (^) 2 



(A^ + n )! 



[a* x a*l 1Sl1 x fa x a] [S2] 



[Ss 



d 7 ™° d 7 ^ (7172) " [1 + 7172(a * ' a)] 

a| [7rf ® tt+J jf (er) 2 |c*) = 

TV + n )(N + n - l)M 2 N no ( j^ ^, [a* x a*]jf 



JV+no-2 



(AT + no)! 



d"° d™° 



d 7 r° d 72 1 

a|(«rt) 2 [ir<S^\a) 



757 (7 2 [1 + 7i 72 (a • a)] ) 



(AH) 



U 2 



N + n )(N + n -l)Af 2 NnoWT ^_ 



[a x a]Jf 



dTTdT 7 / 



(a|(<r + ) 2 (o-) 2 |a) = AT(AT-1) 



;r{7i 2 [l + 7i72(a*-«)] N+ " n - 2 } 
Ml 



J Nn 
Af 2 



(56) 



These equations also give us the mapping of n n and n 2 as 
special cases. For completeness we also give the mapping 
of «4 , which is 



(a| n 2 + 2 V" ttItt^tt 3 + ■k\-k^'k\-k % it 3 n k \a 

= (a| 3 |[ TT* ® TT 1 ]" [ 7T(g> 7t]q 

+ V5 



7T^ (g) 77^1 



(a\ ^2 ■K\^\n\n % n J n k + n x \ol) 



Ot) 



(57) 



Note that the mapping is more complicated than when 
no Pauli exclusion principle is taken into account (no = 
0), due to the distinct property of a m . (It is not a tensor 
anymore.) Also note that 



[«4" = ^sB-d" 



\a x a 



l[0] 



(58) 



Thus, the sum of both is real. Because they always ap- 
pear in a sum in the expectation value of the Hamiltonian 
with respect to the coherent state, the expectation value 
is always real. This is a remarkable sign of consistency. 

The mapping, concerning the individual clusters, is the 
same as given in [l7|. There, one has to take into account 
that the coherent state acquires the form of a direct prod- 
uct of the state describing the relative motion and the one 
giving the cluster coupling 



a) [Ci x C 2 Y 



(59) 



where the last factor refers to the coupling of the two 
cluster states, which is fixed [I?). 

Next we have to expand the above expressions in pow- 
ers of (a ■ a). In [TtJ the no was neglected compared to 
N. Here, we will take into account the contributions of 
no . The list of expansions is 



[1 +7i72(a* • a)] 



N+n 



171=72=1 







( N+ r) 


[(fc-no)lj 



n 2 



{a* -af 



Ana AJio 

d7i"° d 72 "° 

N+n 

£ 

k—riQ 

a "0 a »q 7i72 1 + 7i72 (a -a) 
d 7l d7 2 

iV+nn-l / 

j2 f N + n o - 1 

k—max{nQ — 1,0) 

x (a* ■ a) k 

d"° d"° M JV+n -i 
72 [l + 7x 72 (a -a)J 



(fc + 1) 
(fc + l-no)! 



71— 72 — 1 

I ^ 2 



d7i"° d 72 " 

N+nn-1 / 

fN + n -l 

k—nQ 

x (a* ■ a) k 



(A -no)! 



'71=72=1 

(fc + 1) 
(fc + l-no) 
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d7? 



o(7i72) [1 + 7i72(a* ■ a)} 



N+no-2 



N+n -2 , 
k— max (no —2,0) 

x (a* ■ a) k 



(fc + 2)! 



171=72=1 

n 2 



(fc + 2 - n ) 



H„-(72) 2 [l + 7i 72 («*•<*)] 



d"° cT 
dTT d72 

Af+rio-2 

E 

fc— max(riQ : 0) 

x (a* - a)* 1 



JV+no-2 



'71=72 = 1 



N + no-2 
k 



(fc 



(fc + 2)! 
- ^0) 



fc! 



(fc - n )! 
(60) 



For example, the leading term in the first expression 
is ^^-""^ reo! (a* ■ a) n " which leads in lowest order in 
(a* • a) to the normalization, as defined in Eq. (8) of 
Ref. [131 ■ We will, however assume that N » no >> 1, 
otherwise the resulting expressions are too involved. 

The equations in (|60[) can be simplified, using the ab- 
breviation 



which leads then to Eq. (f3TH) . 
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